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February 2002 campaign at the High Frequency Active Auroral Research Program 
(HAARP) heating facility are analyzed. We find that during injections toward magnetic 
zenith (MZ) the green and red lines gain ~5 R within ~1 s and ~20 R within ~10 s, 
respectively. We term this period the onset of the HF-induced airglow. A model of the 
onset at magnetic zenith is developed. It accounts for background photoelectrons and 
dissociative recombination of 0\. It is shown that heating and acceleration of background 
electrons dominate the airglow onset. We propose a scenario for the generation of strong 
Langmuir turbulence for injections outside the Spitze region, including magnetic 
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1. Introduction 

[2] A distinctive feature of HF modification experiments 
is the excitation of airglow at 630.0 and 557.7 nm by high- 
power, high-frequency (HF) radio waves [e.g., Sipler et al, 
1974; Bernhardt et al, 1989; Pedersen and Carlson, 2001; 
Gustavsson et al, 2001, 2002; Kosch et al, 2000, 2002]. 
Enhancements up to ~500 R (Rayleighs) with the green-to- 
red ratio C, gr as high as >0.3 have been reported 
[Gustavsson et al, 2002; Kosch et al, 2002; Pedersen et 
al, 2003]. The excitation energies e x of the O('D) and 
0(}S) states responsible for the red and green lines are e r = 
1.96 eV and e g = 4.17 eV, respectively. The population of 
energetic, e > Ex, electrons can increase significantly due to 
stochastic and resonant interactions with plasma turbulence 
generated by a heating wave [e.g., Gurevich et al, 1985; 
Dimant et al, 1992; Mantas and Carlson, 1996; Gurevich 
and Milikh, 1997; Istomin and Leyser, 2003]. The former 
raises the electron temperature T e , the latter accelerates a 
group of electrons in the high-energy tail of the initial 
distribution. Both effects are well documented [e.g., 
Carlson et al, 1982; Gustavsson et al, 2001]. 

[3] Gustavsson et al. [2002] emphasized that to interpret 
> 0.1 in terms of electron heating requires unrealizable 
T e > 2 eV, pointing out the importance of electron acceler- 
ation. It is commonly believed that electrons are most 
efficiently accelerated by Langmuir (1) turbulence. The 
generation of Langmuir waves is usually described in terms 
of nonlinear instabilities of ordinary (o) mode pump waves, 
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known as the parametric decay (PDI) or oscillating two- 
stream (OTSI) instabilities [Fejer, 1979]. In the presence of 
background suprathermal (e > T e ) electrons, many more 
energetic electrons are accelerated than would be in Max- 
wellian plasmas [Mishin and Telegin, 1986]. 

[4] Recent observations show that the airglow maximizes 
during HF injections toward magnetic zenith [Kosch et al, 
2002; Pedersen et al, 2003]. The same is true for the 
intensity of Langmuir waves [Isham et al, 1999] and 
electron heating [Rietveld et al, 2003] observed by the 
EISCAT UHF radar. Furthermore, the red line is excited at 
magnetic zenith even at extremely low effective radiative 
power P ~ 2 MW [Pedersen et al, 2003]. 

[5] Ordinary mode waves with incident angles outside 
the Spitze region, 6 > 9 C , reflect at altitudes J7 9 below the 
standard reflection altitude H where the local plasma 
frequency f p ~ 10 4 ,/«7 Hz equals the driver frequency f . 

Here 6 C = strcsin(^J^^sinx), f c is the local electron gyro- 
frequency, n e is the electron density in cm -3 , and x is tne 
magnetic dip angle [e.g., Mjelhus, 1990]. At HAARP \ - 
14.5° and f c ~ 1 .4 MHz at altitudes near 200 km, so that for 
/o = 7 MHz the Spitze angle is 6,. ~ 5.9°. 

[6] Figure 1 shows a schematic of ray trajectories for 
ordinary HF waves injected vertically and toward magnetic 
zenith into a horizontally stratified ionosphere at HAARP. 
Obliquely incident radiation does not form standing-wave 
patterns and swelling is absent. Thus, given P = 150 MW 
and distance R = 250 km, the wave amplitude is E Q ~ 
A.ly/P^IR ~ 0.2 V/m. With/o = 7 MHz and T e = O.ljeV, 



the HF energy density at the reflection point is W = 



io~Vz 



sufficient to drive the PDI/OTSI at H [Fejer, 
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Figure 1. Schematic of ray propagation for ordinary 
HF waves injected toward local vertical and magnetic 
zenith. The magnetic field B direction is indicated by 
arrows. The heights of reflection and upper hybrid 
resonance are shown by horizontal lines. A light dashed 
line shows the Spitze angle direction. The regions of 
excitation of Langmuir turbulence appear as bold lines near 
the reflection points. 

1979]. However, mismatch of frequencies at H & <H suppress 
these instabilities for injections outside the Spitze region. 

[7] Gurevich et al. [2002] suggested that decreased 
plasma densities within striations generated by heating 
waves permit necessary phase matching [cf. Muldrew, 
1978]. Changes in the refraction index (self-focusing) due 
to striations develop within tens of seconds after turn-on and 
explain some features of the spatial distribution of the HF- 
induced airglow. However, the rise time of Langmuir waves 
is ~10 ms [Isham et al, 1999]. Besides, striations are 
generated in the upper hybrid layer [e.g., Vaskov et al, 
1981; Lee and Kuo, 1983]. Hence the reflection height of the 
heating wave must be at or above this height. For magnetic- 
zenith injections at HAARP, this can be satisfied only if 
f < 5.4 MHz. However, the strongest airglow at HAARP 
occurred with higher fo values [Pedersen et al, 2003]. 

[8] Kuo etal. [1997] showed that Langmuir waves can be 
excited by upper hybrid (uh) waves with amplitudes exceed- 
ing the threshold value ~ 0.1 5 V/m. Such w/i-waves can 
be generated through the linear conversion of the o-mode on 
pre-existing field-aligned irregularities [Wong etal, 1981] or 
through the parametric decay o — > uh + Ih [e.g., Istomin and 
Leyser, 1995]; Ih stands for lower hybrid waves. Linear 
conversion within the upper hybrid layer proceeds with no 
set threshold. However, the parametric decay develops ifiso > 
Ef ^ I.6/0 mV/m, provided that the inequality 0.015 < 
\fo — s • f c \ < 0.5 MHz (integer s > 3) is valid. Here / 
stands for the heating frequency in MHz. The rise time of 
the uh-wave is j ull ~ 1-3 ms. Consistent with observations 
[Isham etal, 1999], this makes the generation of Langmuir 
turbulence possible within ~10 ms. 

[9] Besides electron impact, dissociative recombination 
(DR) of O2 is an efficient source of the 0( l D) and 0('5) 
states. Since the rate of dissociative recombination decreases 
with T e , it is usually not considered in the theory of the HF- 
induced airglow. However, the quantum yield for the 0( l S) 



state from dissociative recombination grows significantly 
with T e and with the vibrational temperature of oxygen ions 
[Guberman, 1997; Peverall et al, 2001]. Furthermore, 
charge exchange + + 2 — > 0\~ + O is the dominant source 
of O2 production in the nighttime F region. Its rate increases 
significantly whenever 2 is vibrationally excited [Viggiano 
and Williams, 2001]. We expect a high degree of excitation 
of molecular species in the HF-illuminated region and thus 
enhanced yields of the green line emissions. 

[10] A consistent theory of the HF-induced airglow 
accounting for the ion/neutral chemistry, modification of 
the heated spot, and self-focusing has yet to be worked out. 
We emphasize that responsible processes develop within tens 
of seconds. However, considerable airglow appears within 
first few seconds after the HF transmitter turns on. We 
designate this interval as the onset of HF-induced airglow. 

[11] This paper develops a model for the onset of 
airglow enhancements at magnetic zenith accounting for 
the roles of O2 dissociative recombination and energetic 
photoelectrons. The following section describes the onset 
characteristics from the HF heating experiments at 
HAARP. Section 3 describes the model. In particular, we 
propose a scenario (section 3.3) for the generation of 
strong Langmuir turbulence outside the Spitze region 
including magnetic zenith with upper hybrid waves as 
the primary source. The final section compares modeling 
results with optical measurements. 



2. Airglow Onset: HAARP, February 2002 

[12] The HAARP facility is located outside of Gakona 
Alaska (62.4°N, 145.1 5 °W). During the course of HF 
heating experiments between 03:00 and 05:00 UT in Feb- 
ruary 2002 several passes of Defense Meteorological Sat- 
ellite Program (DMSP) satellites flew to the east and west of 
the HAARP location. Each of these spacecraft carries a pair 
of upward looking particle spectrometers designed to mea- 
sure fluxes of downcoming electrons and ions with energies 
between 30 eV and 30 keV. Consistent with prevailing quiet 
geomagnetic (Kp = 2) conditions, the equatorward boundary 
of auroral precipitation was several degrees in latitude 
poleward of Gakona. However, both before and after the 
local time of sunset, the spectrometers detected fluxes 
of Rdowncoming electrons whose spectra monotonically 
decreased with energies between 30 and 100 eV. These 
subauroral fluxes consist of photoelectrons that originated 
in the still sunlit southern ionosphere. 

[13] Intense green-line emissions were observed between 
03:48 and 05:00 UT (17:48-19:00 LT) on February 13, 
2002 events of the HAARP optics campaign. Readers are 
referred to the report of Pedersen et al. [2003] for graphic 
examples of HAARP- induced airglow and the intensities of 
red/green-line emissions recorded during the period of 
interest. O-mode waves were injected toward magnetic 
zenith at f = 7.8 MHz and at full power of 0.94 MW 
(P — 1 65 MW). The transmitter was programmed to turn on 
for 5 minutes exactly on the minute. Subsequent pulses 
followed 5 minutes pauses. However, the fourth pulse in the 
sequence started at 04:21:18 UT. A transmitter problem 
produced a false start, 04:30:00-04:31:00 UT, prior to the 
long pulse that began at 04:32:12 UT. In the course of this 
long pulse, the critical frequency of the F layer dropped 
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Figure 2. Heater-induced green (top panel) and red 
(bottom panel) lines during February 13, 2002 injections. 
Vertical solid arid dash-dotted lines indicate heater turn-ons 
and turn-offs, respectively. 

below the transmitter frequency and the sequence was 
terminated. 

[14] To account for natural (background) variations, a 
polynomial fit was made to intensities measured along the 
magnetic meridian in each all-sky image, with the heated 
area blocked out. Throughout the studied period the differ- 
ence between the background airglow b\ and that of the 
heated volume 7x> during heater-off periods, was less than 
~3-5 and ~7— 10 R for green- and red-line emissions, 
respectively. This small difference determines the accuracy 
of heater-induced airglow measurements A x = I\ — b\. 

[is] Figure 2 presents a detailed view of the last three HF- 
induced enhancements of February 13, 2002 at f = 
7.8 MHz. Data points are from two red line exposures 
compiled at and 24 s after the beginning of each minute, 
and green-line exposure beginning at 12 s into the minute. 
After 04:00 UT, the exposure duration was 7.5 s. By 
examining cases when the heater turned on part way 
through an exposure, we see that the green-line emissions 
increased by ~5 (~10) R within ~1.5 (~7.5) s, while the 
red line gained ~25 (~30) R within ~13.5 (~20.5) s. 

[16] Starting at 04:27:15 UT on February 9, 2002, a P ~ 
124 MW wave at a frequency 6.8 MHz was injected toward 
magnetic zenith. Figure 3 shows red- and green-line emis- 
sions before and during the pulse turn on from the heated 



volume and the background region. These data show that 
the difference between natural green-line emissions from 
magnetic zenith and ~70 km apart was less than 1 R. The 
same is true for the red-line emission before the heater 
turned on. This suggests that the background airglow was 
nearly uniform. When the heater turned on, green (red)-line 
intensities increased by ~5 (~30) R within ~4.5 (~16.5) s. 
It is relevant to note that Kosch et al. [2002] reported 
~500 R enhancements of the red- and green-line emissions 
produced by the EISCAT superheater (P ~ 600 MW) 
within the first 5-s frame. 

3. Modeling Airglow Onset 

[n] For modeling purposes we assume that the major 
species constituents of the neutral atmosphere and iono- 
spheric plasma are well described by the MSIS-E [Hedin, 
1991] and IRI [Bilitza, 2001] models, respectively. Although 
photodissociation of 2 accounts for ~50% of the back- 
ground airglow, it is unaffected by heating and hence is not 
included. Basic processes and their rate coefficients are 
discussed at length by Solomon et al. [1988], Rees [1989], 
Witasse et al. [1999], and Fox and Sung [2001]. Unless 
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Figure 3. Heater-induced green (top) and red (bottom) 
lines during February 9, 2002 injections. Dots and asterisks 
stand for emissions from the magnetic zenith and back- 
ground regions, respectively. Vertical dashed line indicates 
the heater turn-on time. 
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otherwise noted, we use the rate coefficients of Fox and 
Sung [2001]. 

[is] As mentioned above, radiative emissions from a 
heated volume reflect the superposed results of complex 
wave, plasma, and chemical interactions. To render our 
modeling of this complex chain of interactions intelligible, 
we divide this section into three parts. The first subsection 
concerns the basic contributors to the natural and artificially 
excited radiation budgets, the subsection on the back- 
ground ionosphere estimates the distributions of ion species 
and energetic photoelectrons present at heater altitudes prior 
to turn on. The third subsection describes a flow chart 
encapsulating our concept of how heater-injected energy is 
transformed into the various wave modes that heat and/or 
accelerate ambient electrons that interact with ambient 
neutrals to produce onset green- and red-line emissions. 

3.1. Basic Processes 

[19] The red (r) and green (g) line-photons are emitted by 
atomic oxygen in the transition from the 0('£>) and 0( l S) 
states to 0( 3 P) and 0( l D) states, respectively. In photo- 
chemical equilibrium, the volume emission rate r|x is 
calculated from 



rff =A-[Ox]=A 



Here [O x ] stands for the density of 0('Z>) or 0( l S) in cm -3 ; 
Qx and L\ are the corresponding production and loss rates, 
respectively; A r ~ 0.007, ~ 0.009, A g ~ 1 .2, and A& ~ 
1.3 are the Einstein transition probabilities in s -1 . The 
column emission rate (in R) is found by integrating (1) 
along the line of sight 



T 



■n x (z)& 



(2) 



[20] The rates of electron impact excitation O + e — » 
0( l D, l S) + e and the Schumann-Runge dissociation 2 + 
e — > O + O^D) + e are calculated as follows 



/•OC 

ej = [0]-*£ = [OMW °x(e)*(e)«fe 



(3) 



Here 3>(e) is the electron differential number flux, x\, £x> 
and o\ are the rate coefficients, threshold energies, and 
cross-sections, respectively. We employ electron impact 
cross-sections suggested by Majeed and Strickland [1997]. 
Figure 4 shows three of the major electron impact cross 
sections. It is worth noting that x e r (T e ) and x| (T e ) calculated 
with a Maxwellian distribution of thermal electrons are 
close to the approximations used by Mantas and Carlson 
[1996] and Gurevich and Milikh [1997], respectively. 

[21] The dissociative recombination of 0\ produces 
Q\ R = "x*« e [02]. ^ ne rate coefficients are usually approx- 
imated as %\ R = px* • ot, where (3x* are me quantum yields. 
Estimated yields for 0( l S) were subject to substantial 
discrepancy in the literature for many years. Recent 
experiments [Kella et al., 1997; Peverall et al, 2001; 
A. Petrignani et al., Vibrationally resolved rate coefficients 
and branching in the dissociative recombination of O2, 
submitted to Journal of Chemical Physics, 2003] and theory 




Electron energy (eV) 

Figure 4. Cross-sections of the O^D) and 0( l S) excita- 
tion and the Schumann-Runge dissociation by direct 
electron impact. Squares and dots are the measured cross- 
sections from Doering and Gulcicek [1989]; solid lines are 
interpolations by Majeed and Strickland [1997]. 



(1) [Guberman, 1988, 1997] appear to have reconciled this 
controversy. They established that ct and Px* depend not 



only on T e but also on the vibrational population n^l - 
[OiUWIOi] of 02, where j = 0, 1,. . . designates 2 the 
vibrational level. 
[22] For the ground vibrational state, $ 0) ~ 1.15, ct (0) ~ 2 

10-" 



- 7 t7 63 [cf. Fox and Sung, 2001], and 



of— P (-f (.-!£)) W 

for 1 < r e < 10, where T e = r e [K]/300. Importantly, the 



inequalities ^> > and (3J. 1>2) < # 0) hold. For the 
Maxwell-Boltzmann distribution of N^l with vibrational 
temperatures Tv ~ 0.25 eV, one obtains 2 pf* ~ 1.28^ 0) and 
[3?* ~ 0.85p(°>. 

[23] Collisional quenching # + X — ► O + X is important 
for the 0('D) state, while for 0( l S) at altitudes >150 km it 
can be ignored. Here X stands for O, 2 , N 2 , or e, h . The 
corresponding loss rates are L? = 10~" nx[X], where xS — 
0.65 t° 14 , %S=± 3.2 exp (0.225/t„), x N q 2 ~ 1.8 exp (0.36/t„), 
and x| ~ 28.7 t c 91 . Figure 5 shows the altitude profile of 
the total loss rate calculated with the MSIS-E and IRI 
parameters at 04:30 UT, February 13, 2002. 

[24] Finally, the radiative cascade (RC) 0( l S) -+ OCD) + 
hv 5577 yields Qf c = 1^. 

3.2. Background Ionosphere 

[25] Simultaneous observations from a digisonde located 
at the HAARP site were used to correct the IRI model and 
to determine that the reflection height at magnetic zenith 
decreased from ~250 to ~235 km between 04:00 and 
04:30 UT on February 13, 2002. At the same time, the 
shadow height increased from ~150 to ~230 km indicating 
the presence of photoelectrons [e.g., Doering et al., 1975]. 
Furthermore, the southern-hemisphere region magnetically 
conjugate to HAARP remained in sunlight, even after local 
sunset. Thus a large fraction of the photoelectron spectrum 
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0( 1 D) loss rate, s 1 
Figure 5. 0( l D) loss rate. 



generated in the conjugate ionosphere [e.g., Peterson et al, 
1977] had access to the ionosphere above HAARP. Electron 
spectrometers on DMSP satellites observed the high-energy 
tail of the conjugate photoelectron spectrum. 

[26] Figure 6 shows altitude profiles of the background 
electron temperature T and electron and + densities 
calculated from the corrected IRI model. The twilight 0\ 
density above ~200 km is defined by the charge exchange 
and dissociative recombination 

Here the rate of charge exchange x ce after Viggiano and 
Williams [2001] is 



Figure 7. First, upper hybrid waves are generated due to 
the parametric decay o —* uh + Ih. Its threshold field is 
Eo h ~ 0.08 V/m if the matching conditions are met at 
x u h = kl h rl < 1, where r e is the thermal electron gyrora- 
dius and k uh is the w/i-wave vector. Landau damping of 
short-scale lower hybrid waves raises the threshold other- 
wise [Mishin et al, 1997]. For/ = 7.8 (6.8) MHz and 9 ~ 
14.5°, the matching conditions at the reflection height Z/ 
are easily satisfied for x ull ~ 0.3 (~0.15). Finally, the free 
space field of the incident wave E 0.25 (0.2) V/m 
exceeds E§ h . 

[29] When the amplitude of the primary wA-wave E uh 
exceeds ~10 mV/m, it excites other, lower frequency, 
w/z-wave with the growth rate of order y uh ~ [Zhou 
et al., 1994]. The same is true for subsequently generated 
w/j-waves. The frequency-step of this spect ral tran sfer is 
quite small |6co u /,|/u)„ A ~ 0Ax uh bk uh /k u i, ~ 0.2^/m e /mi, where 
mjnti represents the electron/ion mass ratio. As many 
spectral steps occur before the parametric instability satu- 
rates, the resulting uh-energy spectrum consists of a large 
number AS; 10 of spectral peaks, each of order Eo/(y uh T uh ) 
(well-known weak turbulence "cascading" [e.g., Sagdeev 
et al, 1991]). Thus the total wA-energy density can be 
estimated as W*lh. ~ W ~ 7.5 (5) eV/cm 3 with the r.m.s 
amplitude E* ~ ~ \/AE - 

[30] Upper hybrid waves with amplitudes E uh exceeding 
the threshold value E' uh ~ 0.1 lx\£ V/m generate short- 
scale, k, > 3~ l/2 /r e , Langmuir waves [Kuo et al, 1997]. 
Due to a weak dependence of the growth rate 7/ on k uh , 
one can take E uh ~ is*. As this value greatly exceeds the 
threshold E l uh , the growth rate can be evaluated as 7/ ~ 
v e E uh IEu h [Kuo et al, 1997], which amounts to 7, ~ 
v e y/\Q 2 x u hhEQ- Here v e is the electron elastic (transport) 
collision frequency. 

[31] The growth of Langmuir waves at TJT t < 4 saturates 
via induced scattering by ions that transfers energy toward 



* ce ~ x<») + »^[4.91g(T,„ + 1) - 2.5]) 



(6) 



*»(e) 



•0.5 



4tc s 
3 and p s ~ 



(7) 



3. 



3.3. HF-Perturbed Ionosphere 

[28] Our scenario for the excitation of plasma turbulence 
and subsequent electron energization in the HF-illuminated 
region at magnetic zenith is represented schematically in 



400 



350 



a is the recombination rate in cm /s, n \ J ' = [O 2 0')]/[O 2 ], 
[Y] stands for the density of species Y in cm" , and t ,„ = _ 
1.5t, + t„ - 1 with t,-„ = 7]-„[K]/300. * 
[27] The IRI model predicts T e =T < 0.2 eV (Figure 6). -g 
At these temperatures and n e > 10 5 cm -3 , a Maxwellian .-2 
distribution F M is a good representation of the thermal ~~ 
electron distribution function (F, h ) in the F region [Mishin 
et al, 2000]. The differential number flux ^(e) of photo- 
electrons in the F region [e.g., Rees, 1989] at e > e s = fynvl 
~ 5 eV can be approximated by a power-law function 
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Figure 6. Altitude profiles of the electron temperature 
(solid lines), electron density (dashed line), and + density 
(dotted line). T , T u and T 2 stand for 0, 1, and ~1.5 s after 
turn-on. 
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Figure 7. Block-diagram showing synergy of a model of 
HF-induced airglow at magnetic zenith. Shaded boxes 
indicate the steps to be worked out. 



small k [e.g., Zakharov et al, 1976]. This process is 
governed by the pump u/i-wave power. In particular, if the 
width Aki of the parametrically unstable /-wave spectrum 
exceeds Ak c ~ k,vJlh the Langmuir wave energy W t 
accrues in a state with k — > (Langmuir condensate) 
[Zakharov et al, 1976; Zakharov, 1984]. The spectral width 
can be evaluated as Ah ~ 3(fo/fc) 2 W*/(n e T e r e ). Thus, for 
E > 0.15(7 v e ) 1/6 /\/A < 0.2 V/m the "condensate" 
condition Ak t > k,vjy/ is fulfilled. The dynamics of the 
Langmuir condensate is defined by the modulational insta- 
bility and collapse leading to establishment of strong 
(cavitating) Langmuir turbulence [e.g., Galeev et al., 
1977; Zakharov, 1972, 1984]. 

[32] We emphasize that this applies not only at the 
reflection layer Hq but also well below it, wherever the 
matching conditions for the parametric decay o — > uh + Ih 
are met and local values of x„h < 1. Given the plasma 
density profile shown in Figure 6 and f ~ 6-8 MHz, the 
altitude extent of this region is ~10 km. One should 
compare this value with a few 100-m size of the Airy 



pattern at H , where the PDI/OTSI instabilities develop 
during injections within the Spitze region [Fejer, 1979], We 
believe that this is the key factor in exciting strong airglow 
at magnetic zenith. 
3.3.1. Heating 

[33] Collisional damping of high-frequency o-, uh-, and 
/-waves is the major source of (stochastic) electron heating. 
Mishin et al. [2000] showed that at [N 2 ] > 10 7 /o cm -3 the 
process of 7V" 2 -vibrational excitation dominates the forma- 
tion of the distribution function of ionospheric electrons 
with energies 3.5 > e > e v » — 1-8 eV, which can be 
represented as follows 



F,h (e > e v ,v, 



x exp 



/2Te 



• 5 *.(£..») 



(8) 



Here k, (e) = Vife)lv ee {z) + 0.5z/T e , where v a lv ee is the ratio 
between electron inelastic and electron-electron collision 
frequencies. A total HF wave energy density W HF -C 
lO~ 2 n T e is assumed. At energies e < the distribution is 
close to F M , while at e > 3.5 eV it may slightly deviate from 
F M due to the 0( l D) and OCS) excitation. 

[34] The electron temperature can be evaluated from the 
energy balance 



l.5n e dT e /dt = T e - v e 5 • n e T e - Vw 4 



(9) 



where q e \\ is the (parallel) electron heat flux, T e ~ v e W HF is 
the volume heating rate, and 6 = (v,//v e ) is the coefficient of 
inelastic losses averaged over the total distribution F th . 

[35] Given v e ~ 500 s _1 , T e is fe40 (25) keV/cm 3 /s. Note 
that to match their optical observations, Manias and 
Carlson [1996] used the values of T e between ~50 and 
~125 keV/cm 3 /s at ~260 km for f„ ~ 5 MHz. Gustavsson 
et al. [2001] used T e ~ 60 keV/cm 3 /s for/ ^ 4 MHz to fit 
the electron temperature profile maximum 4000 K at 
~220 km observed by the EISCAT UHF radar. Neither of 
them accounted for the decrease of 8(7^) due to the deviation 
of the thermal electron distribution from Maxwellian. Taking 
that into account, the calculation results of Mantas and 
Carlson [1996] and Gustavsson et al. [2001] can be scaled 
to the Ti a-profiles in Figure 6, pertaining to our case. 
3.3.2. Acceleration 

[36] Resonant Ih- and /-wave-particle interaction accel- 
erates electrons. Musher et al. [1978, 1986] analyzed in 
detail the dynamics of /A-waves excited by an external 
sou rce. At 1 < TJT t < 4, high-frequency waves, f, h > 
f c y/2m e /mj, are dominated by induced scattering by ions. 
The rate of spectral transfer toward smaller frequencies is 
lih ~ (fc 2 /fihWi h /n T e , where W, h is the energy density of 
lower hybrid waves. Equating the growth rate y uh to yth 
yields the energy density at the saturated state Wf h ln T e ~ 

luhMc ~ 10-". 

[3 7] The d ynamics of //z-waves in the low-frequency, f lh < 
fcy/lmjmi, region is dominated by the lower hybrid 
collapse [e.g., Musher et al, 1978]. The threshold energy 
density is quite low Wf h lnT e ~ x^imjm^ifjfa) 2 < 10~ 6 , 
and is surely exceeded in our case. In the course of collapse, 
the longitudinal and transverse dimensions of a //z-cavity, 
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/|l ~ l±y/mj/xihm e , diminish. Finally, the tt-wave energy in 
a collapsing cavity is absorbed by particles via Landau 
damping [Musher et al, 1978, 1986; Sotnikov et al, 1978; 
Shapiro et al, 1993]. The density of accelerated electrons in 
the magnetic field-aligned high-energy tail of the distribu- 
tion function, determined by assuming resonant particles 
carry away all of the energy pumped into collapsing cavities 
[e.g., Shapiro et al, 1993], is of order <10~ 5 « e . This 
//j-contribution appears insufficient to account for the 
observed HF-induced airglow. It is worth noting that 
//i-collapse also leads to transversely accelerated ions that 
absorb the collapsing energy at approximately the same rate 
as electrons [Sotnikov et al, 1978]. 

[38] Strong Langmuir turbulence consists of an ensemble 
of collapsing cavitons that transfer the wave energy toward 
small scales. Small-s cale, k > w,t/v min , waves are absorbed 
by fast v > v mjn 3> yjT e /m e electrons via Landau damping, 
thereby increasing the population of suprathermal electrons. 
The distribution of accelerated electrons F a (e) = ^-$ a (e) in 
weakly magnetized f c <^f p plasmas can be found from the 
kinetic equation [e.g., Galeev et al, 1977] 




£V 10 l£ min 



E< 1 > 
mm 



1<f 



energy, eV 



Figure 8. Differential number fluxes of photoelectrons 
(solid lines) and accelerated electrons (dash-dotted lines) in 
cnrVster'eV -1 for n s = 1 (line 1) and 10 (line 2) cm -3 . 

and the wave energy flux transferred by collapsing cavitons 
and absorbed by accelerated particles 



dF a 
dt 



V- Ov \ mnv J Wl/v k* 



dF a 
dv 



(10) 



(15) 



where w*/(2it) ~'f p and W k are, respectively, the frequency 
and spectral energy density of Langmuir waves. The latter is 
determined from the requirement for energy balance 



dW k d_ 
dt + dk 



dk 
dt 



= T k W k 



(11) 



Here E km is the wave amplitude in a cavity of the scale 
length / min ~ k~ l ~ v m Jw p when collapse is arrested and 
$ is the differential number flux of ionospheric electrons. 

[41] For a Maxwellian electron distribution, one has e^ in 
~ 207^. When photoelectrons are present, $ (e) — ► $s(e) at 
e ;§> T e , and e min is defined by 



Here 



■ 30(n s T e /W,) 2/5 eV 



(16) 



valid for the energy density of Langmuir. turbulence W t » 
n^T^nJmAMishin and Telegin, 1986]. Given n s = 10 cm -3 , 
n e = 6 10* cm -3 , and W, ~ 10 _4 « e ^, from (16) one gets 
£mi n — 1 5 eV. The number of accelerated electrons exceeds 



10. For energies e < £m in , the 



is the rate of Landau damping; k{t) ~ (t c - t)~ 2 is defined *e photoelectron background in the energy range e > e, 
by the collapse law in the absorption region [Pelleiter, 
1982], and t c is the collapse time. 
[39] The steady state solution of (10)— (12) at e > e min = 



by a factor of &=|§ 
distribution is unaffected. 



number fluxes 



Figure 
of photoelectrons 



^eVmmis a power-law function [Pelleiter, 1982; Shapiro 
and Shevchenko, 1984] that yields the differential number 
flux 



electrons (13) for n s = 1 and 10 cm 3 
and W, 



shows 
(7) and 



differential 
accelerated 



io- 4 «x. 



n e = 6 10 cm 



*.(e) = 



Pa 



■0.5 



4TC m 



(13) 



(p a ~ 0.75). One-dimensional (magnetically field-aligned) 
numerical modeling of electron acceleration by strong 
Langmuir turbulence yields $ a (£||) ~ ejf 1 , consistent with 
the ID scaling law for Langmuir collapse [Galeev et al., 
1983; Wang et al, 1997]. A flat distribution of accelerated 
electrons is consistent with the observations of Carlson et 
al. [1982]. 

[40] The minimum energy e min and the density of par- 
ticles in the high-energy tail n a are determined by the 
boundary condition 



$<.(Emin) = $o(e ra in) 



(14) 



4. Discussion 

[42] Figure 9 shows the components of the equilibrium 
volume emission rates Tj>?(l) for the period near 0430 UT. 
To calculate the contributions of photoelectrons (7) and 
accelerated electrons (13), n s = 10 cm and Wi ~ 10 n e T2 
were chosen, yielding e^in — 15 eV for / = 7.8 MHz. As 
expected from (4), the dissociative recombination contribu- 
tion to the red line is reduced much more than that for the 
green line. Apparently, heated electrons dominate the red- 
line emission, while accelerated electrons are the major 
contributor to the green line. Note a factor of about 2 
difference between calculations with F/JT{) and F th {T{) 
consistent with Mishin et al. [2000] (it is ~3.5 for T e = T 2 ). 

[43] In the heated volume the densities of excited oxygen 
atoms and of 0\ ions grow. During the onset period [O('D)] 
and [O2] remain well below the new equilibrium values that 
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Figure 9. Equilibrium volume emission rate components 
of the red (a) and green (b) emissions: Dissociative 
recombination with [0 2 } Q from the background T e = T 
(thin solid line) and heated T e = T 2 (bold solid line) regions; 
phdtoelectrons (thin dashed line) and accelerated electrons 
(bold dashed line); 0( l S) cascading from the background 
(circles) and heated (dots) regions; thermal electrons with 
T e = T (crosses), Ti (filled triangles), and T 2 (open 
triangles); Maxwellian electrons with T e = 7\ (asterisks). 

are of the order of QJL r and K ce [0 2 Vol(T 2 ), respectively. 
Corresponding time-variations of the contributing compo- 
nents to the total volume emission rate during the onset can 
be evaluated assuming that < t, < min[l/(n e a 2 ),l/i r ] — 
15 s and L,. ~ const(f) at 200-300 km. Here C = t - t and 
the heating time t ~ 1 s. As a result, one gets 



W = [oZ}/[ol] -i * ».ta(r ) - a(r 2 )] • C 

This yields 

T! g (0 ~ rf + T)f n e [a(To) - ot(r 2 )] • C 

n.(0 * {rPU + 0.5^ R [L r + „ e a(r 2 )]6[0+]X 



(17) 



(18) 



[44] Using (18) and Figure 9, the variation of the red- and 
green-line emissions during the onset period are calculated 



from (2). The results are shown in Figure 10 together with 
the observed intensities (section 2). There is good apparent 
agreement with the HAARP observations. 

[45] We emphasize that the onset values of the HF- 
induced airglow amount to ^10% of the saturated values. 
As the reflection heights of 7.8 and 6.8 MHz waves at 
magnetic zenith are below the upper hybrid resonance, 
similar to Figure 1 , striations are not generated. Thus self- 
focusing due to striations seems unlikely. This indicates that 
the ionosphere's parameters in the heated volume are 
significantly modified. Apparently, a complete explanation 
of artificially induced airglow involves a complex chain of 
plasma and chemical interactions within the heated volume 
and has yet to be worked out. 

5. Conclusion 

[46] Observations of airglow at 630 nm (red line) and 
557.7 nm (green line) acquired during the February 2002 
campaign at the HAARP heating facility have been ana- 
lyzed. During injections toward magnetic zenith the inten- 
sities of green and red line emissions gain ~5-10 R within 
~1 s and ~20 R within ~10 s, respectively. We call this 
period the onset of HF-induced airglow and develop an 
explanatory model for the growth of emissions from mag- 
netic zenith. The model accounts for the roles played by 
ambient photoelectrons and dissociative recombination of 
O* ions and shows that heating and acceleration of back- 
ground electrons, possibly by strong Langmuir turbulence, 
dominate the airglow onset. 

[47] We also propose a new scenario for the generation of 
strong Langmuir turbulence for injections outside the Spitze 
region, including magnetic zenith, with upper hybrid waves 
as the primary source of Langmuir waves. Those waves 
accumulate in the region of small wavevectors due to 
induced scattering by ions (Langmuir condensate). Waves 
in the condensate are subject to the modulational instability 
and collapse thus leading to strong turbulence. This scenario 



30 
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20 
rr 15 
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10 
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Figure 10. Calculated intensities of the HF-induced red 
(solid line) and green (dashed line) emissions during the 
onset period. The observed intensities are shown by open 
circles (red) and triangles (green), respectively. 
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applies not only at the reflection layer but also well below it, 
wherever the matching conditions for the parametric decay 
o — > uh + Ih are met and local values of x uh < 1. For 
injections toward magnetic zenith and heating frequencies 
>5.4 MHz, this occurs below the altitude of the upper 
hybrid resonance. The large altitude extent of the strong 
turbulence region appears to be the key factor in exciting 
strong airglow at magnetic zenith. 

[48] Acknowledgments. This research was supported in part by the 
HF Active Auroral Research Program (HAARP) under AFRL contract 
F19628-02-C-0012 with Boston College and by Air Force Office of 
Scientific Research under tasks 23 11 AS and 23 USD. 

[49] Arthur Richmond thanks Bjorn Gustavsson and Michael 
T. Rietveld for their assistance in evaluating this paper. 

References 

Bernhardt, P., C. Tepley, and L. Duncan (1989), Airglow enhancements 
associated with plasma cavities formed during ionospheric heating 
experiments, J. Geophys. Res., 94, 9071. 

Bilitza, D. (2001), International Reference Ionosphere 2000, Radio Sci., 36, 
261. 

Carlson, H., V. Wickwar, and G. Mantas (1982), Observations of fluxes of 
suprathermal electrons accelerated by HF excited Langmuir instabilities, 
J. Atmos. Terr. Phys., 12, 1089. 

Dimant, Y., A. Gurevich, and K. Zybin (1992), Acceleration of electrons 
under the action of intense radio-waves near electron cyclotron harmo- 
nics, J. Atmos. Terr. Phys., 54, 425. 

Doering, J., and E. Gulcicek (1989), Absolute differential and integral 
electron excitation cross sections for atomic oxygen: The 3 P — » '/) and 
3 P -» 'S transitions from 4.0 to 30 eV, J. Geophys. Res., 94, 1541. 

Doering, J., W. Peterson, C. Bostrom, and J. Armstrong (1975), Measure- 
ment of low-energy electrons in the day airglow and day side auroral 
zone from Atmosphere Explorer C, J. Geophys. Res., 80, 3934. 

Fejer, J. (1979), Ionospheric modification and parametric instabilities, Rev. 
Geophys., 17,135. 

Fox, J., and K. Sung (2001), Solar activity variations of the Venus thermo- 
sphere/ionosphere, J. Geophys. Res., 106, 21,305. 

Galeev, A., R. Sagdeev, V. Shapiro, and V. Shevchenko (1977), Langmuir 
turbulence and dissipation of high-frequency energy, Sov. Phys. JETP, 46, 
711. 

Galeev, A., R. Sagdeev, V. Shapiro, and V. Shevehenko (1983), Beam 
plasma discharge and suprathermal electron tails, in Active Experiments 
in Space (Alpbach, Austria), Eur. Space Agency Spec. Publ., ESA-SP- 
195, p. 151. 

Guberman, S. (1988), The production of O('D) from dissociative recombi- 
nation of Ot Planet. Space Sci., 36,47. 

Guberman, S. (1997), Mechanism for the green glow of the upper atmo- 
sphere, Science, 278, 1276. 

Gurevich, A., and G. Milikh (1997), Artificial airglow due to modifications 
of the ionosphere by powerful radio waves, J. Geophys. Res., 102, 1389. 

Gurevich, A., Y. Dimant, G. Milikh, and V. Vaskov (1985), Multiple accel- 
eration of electrons in the regions high-power radio-wave reflection in the 
ionosphere, J. Atmos. Terr. Phys., 47, 1057. 

Gurevich, A., K. Zybin, H. Carlson, and T. Pedersen (2002), Magnetic 
zenith effect in ionospheric modifications, Phys. Lett. A, 505(5), 264. 

Gustavsson, B., et al. (2001), First tomographic estimate of volume distri- 
bution of HF-pump enhanced airglow emission, J. Geophys. Res., 106, 
29,105. 

Gustavsson, B., B. Brandstrom, A. Steen, T. Sergienko, T. Leyser, 
M. Rietveld, T. Aso, and M. Erjii (2002}, Nearly simultaneous images 
of HF-pump enhanced airglow at 6300 A and 5577 A, Geophys. Res. 
Lett., 29(24), 2220, doi:10.1029/2002GL015350. 

Hedin, A. (1991), Extension of the MSIS thermospheric model into the 
middle and lower atmosphere, J. Geophys. Res., 96, 1159. 

Isham, B., M. Rietveld, T. Hagfors, C. LaHoz, E. Mishin, W. Kofman, 
T. Leyser, and A. van Eyken (1999), Aspect angle dependence of HF 
enhanced incoherent backscatter, Adv. Space Res., 24(%), 1003. 

Istomin, Y, and T. Leyser (1995), Parametric decay of an electromagnetic 
wave near electron cyclotron harmonics, Phys. Plasmas, 2, 2084. 

Istomin, Y., and T. Leyser (2003), Electron acceleration by cylindrical 
upper hybrid oscillations trapped in density irregularities in the iono- 
sphere, Phys. Plasmas, 10, 2962. 

Kella, D., L. Vejby-Christenson, P. Jonson, H. Pedersen, and L. Andersen 
(1997), The source of green line emission determined from a heavy-ion 
storage ring experiment, Science, 276, 1530. 



Kosch, M., M. Rietveld, T. Hagfors, and T. Leyser (2000), High-latitude 
HF-induced airglow displaced equatonvard of the pump beam, Geophys. 
Res. Lett., 27, 2817. 

Kosch, M., B. Gustavsson, and M. Rietveld (2002), An overview of high- 
latitude induced aurora from EISCAT, paper presented at World Space 
Congress, Am. Inst, of Aeron. and Astron., Houston, Tex. 

Kuo, S., M. Lee, and P. Kossey (1997), Excitation of oscillating two stream 
instability by upper hybrid pump in ionospheric heating experiments at 
Troms0, Geophys. Res. Lett, 24, 2969. 

Lee, M., and S. Kuo (1983), Excitation of upper hybrid waves by a thermal 
parametric instability, J. Plasma Phys., 30, 30,463. 

Majeed, T, and D. J. Strickland (1997), New survey of electron impact 
cross sections for photoelectron and auroral electron energy loss calcula- 
tions, J. Phys. Chem. Ref. Data, 26, 335. 

Mantas, G., and H. Carlson (1996), Reinterpretation of the 6300-A airglow 
enhancements observed in the ionosphere heating experiments based on 
analysis of Plateville, Colorado, data, J. Geophys. Res., 101, 195. 

Mishin, E., and V. Telegin (1986), Spectrum of suprathermal electrons in 
the auroral plasma, Sov. J. Plasma Phys., 12, 509. 

Mishin, E., T. Hagfors, and W. Kofman (1997), On origin of outshifted 
plasma lines during HF modification experiments, J. Geophys. Res., 102, 
27,265. 

Mishin, E., H. Carlson, and T. Hagfors (2000), On the electron distribution 
function in the F region and airglow enhancements during HF modifica- 
tion experiments, Geophys. Res. Lett, 27, 2857. 

Mjolhus, E. (1990), On linear conversion in a magnetized plasma, Radio 
Sci., 25, 1321. 

Muldrew, D. (1 978), The role of field-aligned ionization irregularities in the 
generation of the HF-induced plasma line at Arecibo, J. Geophys. Res., 
83, 2552. 

Musher, S., A. Rubenchik, and B. Sturman (1978), Collective effects asso- 
ciated with lower hybrid heating of plasma, Plasma Phys., 20, 1131. 

Musher, S., A. Rubenchik, and I. Shapiro (1986), Nonlinear effects in the 
propagation of an ion beam across a magnetic field, Sov. Phys. JETP, 63, 
519. 

Pedersen, T., and H. Carlson (2001), First observations of HF heater-pro- 
duced airglow at the High Frequency Active Auroral Research Program 
facility: Thermal excitation and spatial structuring, Radio Sci., 36, 1013. 

Pedersen, T., M. McCarrick, E. Gerken, C. Selcher, D. Sentman, 
H. Carlson, and A. Gurevich (2003), Magnetic zenith enhancement of 
HF radio-induced airglow production at HAARP, Geophys. Res. Lett., 
30(4), 1169, doi:10.1029/2002GL016096. 

Pelleiter, G. (1982), Generation of a high-energy electron tail by strong 
Langmuir turbulence in a plasma, Phys. Rev. Lett, 49, 782. 

Peterson, W., J. Doering, T. Potemra, R. McEntire and C. Bostrom (1977), 
Conjugate photoelectron fluxes observed on Atmosphere Explorer C, 
Geophys. Res. Lett., 4, 109. 

Peverall, R., et al. (2001), Dissociative recombination and excitation of 0£: 
Cross sections, product yields and implications for studies of ionospheric 
airglows, J. Chem.. Phys., 114, 6679. 

Rees, M. (1989), Physics and Chemistry of the Upper Atmosphere, Cam- 
bridge Univ. Press, New York. 

Rietveld, M., M. Kosch, N. Blagoveshchenskaya, V. Kornienko, T. Leyser, 
and T. Yeoman (2003), Ionospheric electron hearing, optical emissions and 
striations induced by powerful HF radio waves at high latitudes: Aspect 
angle dependence, / Geophys. Res., 108(A4), 1141, doi:10.1029/ 
2002JA009543. 

Sagdeev, R., V. Shapiro, and V. Shevchenko (1991), Anomalous dissipation 
of strong electromagnetic waves in a laser plasma, in Physics of Laser 
Plasma, vol. 3, edited by A. Rubenchik and S. Witkowski, p. 271, North- 
Holland, New York. ■ ■ . 

Shapiro, V., and V. Shevchenko (1984), Strong turbulence of plasma oscil- 
lations, in Basic Plasma Physics, edited by A. Galeev and R. Sudan, 
vol. 2, p. 123, North-Holland, New York. 

Shapiro, V., V. Shevchenko, A. Sharma, K. Papadopoulos, R. Sagdeev, and 
V. Lebedev (1993), Lower-hybrid turbulence at cometary bow wave and 
acceleration of cometary protons, J. Geophys. Res., o 98, 1325. 

Sipler, D., E. Enemark, and M. Biondi (1974), 6300-A intensity variations 
produced by the Arecibo ionospheric modification experiment, J. Geo- 
phys. Res., 79, 4276. 

Solomon, S,, P. Hays, and V. Abreu (1988), The auroral 6300 A emission: 
Observation and modeling, J. Geophys. Res., 93, 9867. 

Sotnikov, V., V. Shapiro, and V. Shevchenko (1978), Macroscopic effects of 
collapse at lower hybrid resonanse, Sov. J. Plasma Phys., 4, 250. 

Vaskov, V. V., A. Gurevich, and A. N. Karashtin (1981), Thermal self- 
focusing instability of plasma waves near resonance, Geomagn. Aeron., 
Engl. Transl., 21, 724. 

Viggiano, A., and S. Williams (2001), Ion-molecule kinetics at high tem- 
peratures (300-1800 K): Derivation of internal energy dependencies, 
Adv. Gas-Phase Ion Chem., 4, 85. 



9 of 10 



A02305 



MISHIN ET AL.: HF-INDUCED AIRGLOW ONSET AT MAGNETIC ZENITH 



A02305 



Wang, J., D. Newman, and M. Goldman (1997), Vlasov simulations of 
electron heating by Langmuir turbulence near the critical altitude in the 
radiation-modified ionosphere, J. Atmos. Sol. Terr. Phys., 59, 2461. 

Witasse, 0., J. Lilensten, and C. Lathuillere (1999), Modeling of the 01 
630.0 and 557.7 nm thermospheric dayglow during EISCAT-WINDII 
coordinated measurements, J. Geophys. Res., 104, 24,639. 

Wong, A., G. J. Morales, D. Eggleston, J. Santoru, and R. Behnke (1981), 
Rapid conversion of electromagnetic waves to electrostatic waves in the 
ionosphere, Phys. Rev. Lett., 47, 1340. 

Zakharov, V. (1972), Collapse of Langmuir waves, Sov. Phys. JETP, 35, 
908. 

Zakharov, V. (1984), Collapse and self-focusing of Langmuir waves, in 
Basic Plasma Physics, edited by A. Galeev and R. Sudan, vol. 2, p. 81, 
North-Holland, New York. 



Zakharov, V., S. Musher, and A. Rubenchik (1976), Weak Langmuir turbu- 
lence of an isothermal plasma, Sov. Phys. JETP, 42, 80. 

Zhou, H., J. Huang, and S. Kuo (1994), Cascading of the upper hybrid/ 
electron Bernstein wave in ionospheric heating experiments, Phys. 
Plasmas, 1, 3044. 



W. J. Burke and T. Pedersen, Air Force Research Laboratory, Space 
Vehicles Directorate, Hanscom Air Force Base, MA 01731, USA. 
(william.burke@hanscom.af.mil; todd.pedersen@hanscom.af.mil) 

E. V. Mishin, Boston College, Institute for Scientific Research, 402 St. 
Clements Hall, 140 Commonwealth Avenue, Chestnut Hill, MA 02467- 
3862, USA. (evgenii.mishin@hanscom.af.mil) 



10 of 10 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the 
data needed, and completing and reviewing this collection of information. Send comments regarding this burden estimate or any other aspect of this collection of Information, including suggestions for reducing 
this burden to Department of Defense, Washington Headquarters Services, Directorate for Information Operations and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202- 
4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comply with a collection of Information if it does not display a currently 
valid OMB control number. PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 


1. REPORT DATE (DD-MM-YYYY) 

13-06-2005 RE PR TNT 




4. TITLE AND SUBTITLE 

On the Onset of HF-induced Airglow at HAARP 


5a. CONTRACT NUMBER 
5b. GRANT NUMBER 

5c. PROGRAM ELEMENT NUMBER 

61102F 


6. AUTHOR(S) 

E.V. Mishin*, W.J. Burke and T. Pedersen 


5d. PROJECT NUMBER 

2311 

5e. TASK NUMBER 

SD 

5f. WORK UNIT NUMBER 

A3 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Air Force Research Laboratory/VSBXP 

29 Randolph Road 

Hanscom AFB MA 01731-3010 


8. PERFORMING ORGANIZATION REPORT 
NUMBER 

AFRL-VS-HA-TR-2005-1061 


9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 


10. SPONSOR/MONITOR'S ACRONYM(S) 

11. SPONSOR/MONITOR'S REPORT 
NUMBER(S) 



12. DISTRIBUTION / AVAILABILITY STATEMENT 

Approved for Public Release; Distribution Unlimited. 



*Boston College, Institute for Scientific Research, Chestnut Hill, MA 



13. SUPPLEMENTARY NOTES 

REPRINTED FROM: JOURNAL OF GEOPHYSICAL RESEARCH, Vol 109, A02305, 
doi: 10.1029/2003JA010295, 2004. 

14. ABSTRACT 

[i] Observations of airglow at 630 nm (red line) and 557.7 nm (green line) during the 
February 2002 campaign at the High Frequency Active Auroral Research Program 
(HAARP) heating facility are analyzed. We find that during injections toward magnetic 
zenith (MZ) the green and red lines gain ~5 R within ~1 s and ~20 R within ~10 s, 
respectively. We term this period the onset of the HF-induced airglow. A model of the 
onset at magnetic zenith is developed. It accounts for background photoelectrons and 
dissociative recombination of 0\. It is shown that heating and acceleration of background 
electrons dominate the airglow onset. We propose a scenario for the generation of strong 
Langmuir turbulence for injections outside the Spitze region, including magnetic 
zenith. 
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